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Neutrinoless double-/! decay {Qvfdf} decay) is a hypothetical process that can occur if the neutrino 
is its own antiparticle. The COBRA collaboration operates a demonstrator to search for these 
decays at the Laboratori Nazionali del Gran Sasso in Italy using CdZnTe semiconductor detectors. 

The exposure of 234.7kgd considered in this analysis was collected between September 2011 and 
February 2015. The analysis focuses on the decay of the nuclides ^^®Te, ^°Zn, ^®°Te and 

^^®Cd. A Bayesian analysis is performed to estimate the signal strength of Qvfljd decay. No signal 
is observed for any of these nuclides. Therefore, the following half-life limits at 90 % credibility 
are set: > 1.6 x lO^^r (^^"‘Cd), > 1.9 x lO^Vr (^^*Te), > 6.8 x 10^® yr (^°Zn), 

T °/2 > 6.1 X lO^Vr (^"^“Te), and > 1.1 x lO^^ yr (“®Cd). 


I. INTRODUCTION 

An open question in neutrino physics is the absolute 
mass and nature of the neutrino. Is it a Majorana or 
a Dirac particle, i.e. is it its own anti-particle or not? 
These question can potentially be answered if neutrino¬ 
less double-/3 decay (0;//3/3 decay) is observed. 

Double-/! decay is a rare second-order weak process 
that changes the atomic number by 2 units, while leav¬ 
ing the atomic mass constant It comes in two modes: 
neutrino-accompanied double-/! decay {2vf3j3 decay) re¬ 
sults in two charged leptons and the corresponding neu¬ 
trinos, (A, Z) {A, Z -b 2) -h 2e“ -I- 2D. This decay mode 
has been observed in several isotopes, and half-lives be¬ 
tween 7.1 X 10^® yr (^°®Mo) and 7.7 x 10^^ yr (^^®Te) have 
been measured |5H1]. Neutrinoless double-/! decay [S], 
which is mediated via an exchange of a virtual neutrino, 
can only occur if the neutrino is a Majorana particle [1]. 
This process is forbidden in the standard model of par¬ 
ticle physics because it violates the lepton number by 2 
units: (A, Z) —>■ (A, Z±2)-|-2e^. The half-life of this de¬ 
cay is inversely proportional to the square of the effective 
Majorana neutrino mass. 

QvfiP decay has been searched for in several isotopes, 
but has not been observed, yet. Recent experiments 
report upper limits on the half-life. The CERDA 
collaboration searches for OvPP decay of ^®Ge using 
high purity germanium detectors; the currently best 
limit on the half-life is > 1.9 x lO^^yr at 90% 

credibility j^. The KamLAND-Zen and EXO-200 
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collaborations have searched for the Oi//!/! decay of 
^®®xe mm- Both experiments have observed zero events 
in their region of interest, which results in a half-life 
lower limit of > 2.6 x 10^® yr for the KamLAND- 
Zen Collaboration and > lA x 10^® yr for the 

EXO-200 Collaboration, both at a 90% confidence 
level. The CUORE collaboration set a new limit of 
'^1/2 > 4 X 10^"^ yr at 90% credibility for Ovfiff decay of 
^®°Te |TD]. The Solotvina experiment sets a limit of 
T //2 > 1-7 X 10^® yr for ^^®Cd [H]. Other experiments 
have searched for the 0vj3f3 decay of ^^^Cd and ™Zn by 
using CdW02 and ZnW02 scintillator crystals. Both 
experiments have observed zero events which results 
in half-life limits of T/)/ > 1-1 x 10^^ yr for ^^“^Cd and 

T //2 > 3.2 X 10^® yr for ^®Zn, both limits are reported 
at 90% confidence level min]. For an overview of the 
current experimental status, see e.g. Ref. |14| . 

This paper reports on the search for Qvfifi decay 
of several isotopes using the COBRA demonstrator. 
The CdZnTe semiconductor crystals deployed in the 
setup contain nine double-/! isotopes with several 
decay modes. A focus is placed on the five /!“/!“ 
ground-state-to-ground-state transitions of the isotopes 
ii4cd, i28Te, ^°Zn, ^^OTe and “^Cd with Q-values 
between 542 keV and 2 813keV. Section II introduces 
the COBRA demonstrator. The data taking and data 
quality criteria are discussed in Section III. In Section IV, 
the energy reconstruction and the calibration procedure 
are explained. The pulse shape analysis is discussed in 
Section V. The event selection for the analysis is dis¬ 
cussed in Section VI. Section VII explains the methods 
used for the signal estimation. Section VIII introduces 
all systematic uncertainties which are considered in the 
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signal estimation. The results are discussed in Section IX 
and Section X concludes the article. 


II. THE COBRA DEMONSTRATOR 

The COBRA demonstrator is located at the Labo- 
ratori Nazionali del Gran Sasso (LNGS) of the INFN 
in L’Aquila, Italy. The setup is described in detail in 
Ref. jT5]. The demonstrator comprises 64 CdZnTe 
coplanar-grid semiconductor detectors arranged in four 
layers of 4 x 4 detectors. Each detector has a size of 
1x1x1 cm^ and a mass of 5.9 g. The detectors are sup¬ 
ported by a polyoxymethylene frame installed in a sup¬ 
port structure made of electroformed copper. The setup 
is constantly flushed with evaporated nitrogen to sup¬ 
press radon-induced background. The detectors are sur¬ 
rounded by 5 cm of electroformed copper, followed by a 
5-cm layer of ultra-low activity lead (<3 Bq/ kg of ^^°Pb) 
and 15 cm of standard lead. A 5-cm borated polyethy¬ 
lene layer is used to reduce the neutron flux inside the 
setup. A charge-sensitive preamplifier integrates the cur¬ 
rent pulses and converts the single-ended detector pulses 
into differential signals in order to suppress electronic 
noise during transmission. After linear amplification, the 
pulse shapes are digitized using 100 MHz flash analog-to- 
digital converters (FADGs) and written to disk. 


III. DATA TAKING AND DATA QUALITY 
CRITERIA 

The data taking of the GOBRA demonstrator started 
in September 2011 after the commissioning of a partial 
setup that contained only one detector layer. The other 
three layers were installed successively in 2012 and 2013. 
The demonstrator was completed in November 2013 with 
64 installed detectors. Gurrently, 61 out of 64 detectors 
operate under stable conditions. For further information 
see Ref. [TB] . 

This publication uses data corresponding to an expo¬ 
sure of 234.7kgd, recorded between September 2011 and 
February 2015. A set of quality criteria based on the 
observed noise level of single detectors is applied. A 
small fraction of data taking runs is discarded because 
of a failure of the nitrogen flushing. In such cases radon 
diffused into the setup and the overall background rate 
increased for the duration of these runs. After the data- 
quality selection, the dataset corresponds to an exposure 
of 216.1kg d. 

For the analysis of ^^"^Gd, runs with energy thresholds 
larger than 421 keV are discarded due to the lowered trig¬ 
ger efficiencies. Therefore, the corresponding exposure is 
reduced to 212.8 kg d. 


IV. ENERGY CALIBRATION AND 
RESOLUTION 

The energy deposited inside a GdZnTe detector is es¬ 
timated by the pulse heights of signals measured with 
electrodes on the detector surface. Because the mobility¬ 
lifetime product for holes in GdZnTe is 2 orders of mag¬ 
nitude lower than for electrons, the hole contribution to 
the induced signal depends strongly on the interaction 
depth. Therefore, GdZnTe detectors are typically de¬ 
signed to have a coplanar-grid configuration, which is a 
modification of the Frisch grid principle used for semicon¬ 
ductors m- This approach is based on two interleaved 
comb-shaped anodes, a collecting anode (GA) on ground 
potential and a non-collecting anode (NGA) on a slightly 
negative potential. The deposited energy is then esti¬ 
mated to be 

E cx Hca — wFInca, (1) 

where Hca and Hnca are the pulse heights of the two 
anode signals and w is a weighting factor introduced to 
compensate for electron trapping. This weighting factor 
is determined and optimized for each detector by cali¬ 
bration measurements. In this configuration, the signal 
shape is only generated by the electron contribution. Be¬ 
sides the improved energy resolution, the coplanar-grid 
technique allows for the reconstruction of the interaction 
depth. An analytical model for the calculation of the 
interaction depth is discussed, e.g., in Ref. |18| . 



Energy [keV] 


FIG. 1: The effective FWHM as a function of the 
energy determined with calibration measurements 
(markers). The used line energies are 238 keV (^^^Pb), 
511 keV (annihilation 7 radiation), 1275keV (^^Na) and 
2 614keV (^°®T1). Also shown is the parametrization 
(red line) including the corresponding uncertainty 
(yellow band). 

A ^^Na and a 22877 ^ source are used to determine the 
energy calibration of each detector. The energy threshold 
of each detector is set depending on the noise level. The 
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FIG. 2: Difference between the estimated and the true 
energy for several 7 lines (markers) as well as the 
systematic uncertainty associated with the energy 
calibration (red dashed line) 


thresholds range from 40keV up to 550 keV for all ROIs 
except the ^^^Cd ROI, where the highest threshold is 
421 keV. Higher noise levels are observed after periods of 
detector maintenance if the detectors have been exposed 
to light or humidity. 

The effective energy resolution of the demonstrator 
setup is estimated from a combined spectrum of all detec¬ 
tors, each weighted with its corresponding exposure. The 
7 peaks from measurements with the calibration sources 
are modeled by double-Gaussian functions from which 
the full width at half maximum (FWHM) is calculated. 
The parametrization of the energy resolution as a func¬ 
tion of the energy takes into account the charge carrier 
production and a noise component. The FWHM takes 
the form 


FWHM(F;) = y/oETb, (2) 


where a is a Poisson component and b is associated with 
the noise. Figure 1 shows the effective FWHM as a func¬ 
tion of the incident energy. It ranges from about 15keV 
at an energy of 238 keV up to approximately 30keV at 


2.6 MeV. 

A linear function is used to determine the energy cal¬ 
ibration of each detector. [Figure 2\ shows the difference 
between the peak position and the calculated calibration 
function of several 7 lines. 


V. PULSE SHAPE ANALYSIS 

The analysis of the pulse shapes is a powerful tool to 
identify energy deposits on the detector surface caused, 
e.g., by a and 13 radiation. Three characteristic quan¬ 
tities are defined in the following, which are used to 
suppress events which deposit energies close to the elec- 
trodes(depth requirement) and lateral surface events 


(LSE requirement). The pulse shape analysis is per¬ 
formed using MAnTiGORE m- The recorded pulse 
shapes can be divided into four time intervals: The pre¬ 
pulse baseline, the charge-cloud drift, the charge collec¬ 
tion, and the postpulse baseline. The prepulse and post¬ 
pulse baselines of the difference pulse are used to deter¬ 
mine the pulse height of the signal, which is proportional 
to the deposited energy. During the charge-cloud drift 
the induced signal is approximately the same on both 
anodes. The difference signal is therefore approximately 
zero. Once the charge cloud approaches the anode side 
the grid bias potential becomes important. The electrons 
are collected by the GA, leading to a sharp rise of the GA 
signal and a sharp fall of the NGA signal. The resulting 
weighted difference signal is a function of the energy and 
does not depend on the interaction depth. The interac¬ 
tion depth is reconstructed by the ratio of the cathode 
pulse height and the difference pulse height. The cathode 
signal itself is reconstructed from the sum of both anode 
signals. 

If the energy is deposited on the lateral surfaces near 
the GA, the difference signal shows an earlier rise than 
the bulk signal. This behavior is characterized by the 
early rise time (ERT) value, which is defined as the time 
difference between the 3% point to 50% point of the 
maximum pulse height. On the NGA near lateral sur¬ 
faces, the difference pulse shows a small dip below the 
baseline before the sharp rise of the signal. It is charac¬ 
terized by the dip value, which is defined as the difference 
of the baseline and the pulse minimum in a window of 
30 samples ending at the 50 % pulse height point on the 
right edge. More information about these quantities are 
given in Ref. |20| . 


VI. EVENT SELECTION 

Pulse shapes are recorded if at least one detector 
records an energy deposition above its energy threshold. 
Before any software requirement is applied, the spectrum 
between 400 keV and 5 MeV contains 166 668 counts; 
104 013 counts remain after requiring the reconstructed 
interaction depth to be in the interval [0.2,0.97] and thus 
removing events close to the cathode or the anode side. 
Requiring the dip and ERT value below 53 and 8 , respec¬ 
tively, reduces the number of events to 65 176. For the 
current analysis, no anticoincidence between the detec¬ 
tors is required. This will be studied in the future. 

[Figure 3] shows the spectrum between 400 keV and 
5 MeV after the event selection. Three 7 lines are iden¬ 
tified. The e+e” annihilation line at 511 keV is most 
likely dominated by ^^Na, which is also responsible for 
the line visible at 1 275 keV. A primordial source of radi¬ 
ation is causing a line at 1460 keV. The small peak at 
3.5 MeV is most likely caused by a-radiation on the cath¬ 
ode side. They are induced by Radon remnants which are 
attracted to the cathode after short failures of the nitro¬ 
gen hushing. The detectors of the fourth layer have a 
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coated cathode side with a thickness of 30 ptm. This is 
not sufficient to stop all a particles, but instead reduces 
their kinetic energy. 



FIG. 3: The energy spectrum of all detectors after the 
full event selection. Three 7 lines are identified: the 
annihilation line at 511-keV, the 1 275-keV line of ^^Na 
and the 1460-keV line of The line at 3.5 MeV is 
most likely caused by a radiation on the cathode side. 


VII. BACKGROUND MODELING AND SIGNAL 
ESTIMATION 


Lines from 7 radiation and the signal process have in¬ 
trinsic widths smaller than the energy resolution. The 
energy resolution varies from detector to detector. Thus, 
a 7 line in the dataset combining all detectors should in 
principle be described by a superposition of 64 Gaussian 
peaks. For simplicity such lines are described by double- 
Gaussian functions which model the energy resolution of 
the full detector array: 


^line / s 
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(5) 

where is the expected number of counts for an 

energy interval Aif^, A is the amplitude, F is the fraction 
of the first Gaussian, Eq is the peak position and cti 2 are 
the widths of the Gaussian functions. A can be the total 
number of signal events, S, or Siine the total number of 
line-shaped background events. Studies show that cti = 
2(72 is a reasonable choice which describes the full energy 
peaks well and ensures that the fraction of both Gaussian 
functions is stable over the whole energy range. 

The parameters i3. A, S and, if applicable, i?iine are 
free parameters of the model, for which uniform prior 
probabilities are assumed. The parameters F, (Ti 2 and 
Eq are nuisance parameters, which are used to model the 
corresponding systematic uncertainties, assuming Gaus¬ 
sian prior probabilities. 

The half-life of Qvpfi decay is calculated from the signal 
strength as 


A Bayesian analysis m is performed to estimate the 
signal strength of potential Qvfifi decay. For each isotope 
a ROI is defined as as the ±3FWHM interval around 
the Q-value. For ^^^Cd and ^^°Te the ROI is extended 
to cover the 7 lines which can contribute to the back¬ 
ground. The statistical model is a binned likelihood in 
the energy distribution describing Poisson fluctuations of 
the number of events ni in each bin, i.e. 


n 

p(ri|P,,f4,J7ine) = 

2 = 1 


o-i^l+^l + ^hne) 




(3) 

with being the expectation of the signal, the expec¬ 
tation of the continuous background in each bin AEi and 
Wine expectation of a line-shaped background, which 
is only greater than zero if a line is present in the ROI. 

Based on previous GOBRA data, the background is 
modeled by an exponential function: 
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where B is the total number of background events in the 
ROI, Rmin/max the interval boundaries of the ROI 
and A is the decay parameter of the exponential function. 
The comparison with a background model using a second 
order polynomial results in a difference in the estimated 
half-life of less than 2 % 


rrOvpp _ ln( 2 ) NXe 
^ 1/2 — - ’ 

with N being the number of atoms of the evaluated iso¬ 
tope per kilogram detector mass, X is the total exposure, 
and s is the observed signal strength of the decay. 

The total efficiency e includes the intrinsic efficiency for 
the signal process and the efficiencies of the pulse shape 
analysis requirements. The systematic uncertainties as¬ 
sociated with the efficiency e and the number of atoms N 
are also modeled using nuisance parameters and assum¬ 
ing Gaussian prior probabilities. 


VIII. SYSTEMATIG UNCERTAINTIES 


The systematic uncertainties discussed in the following 
are used as the standard deviations of the Gaussian priors 
for the corresponding nuisance parameters. These are 
summarized in ITableTl 

The uncertainties of the Q-values range from 0.01 keV 
to 2.1 keV [55]. The systematic uncertainty associated 
with the energy calibration, and thus on the peak posi¬ 
tion Eq, is estimated by the largest difference between 
the calibration function and the means of Gaussian fits 
of all 7 lines. The maximum difference is 1.9 keV. 
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The total uncertainty associated with the mean of the 
0v/3(3 decay line is calculated using uncertainty propaga¬ 
tion of Q-value uncertainty and the uncertainty of the 
energy calibration. 

The uncertainties of F and ai are estimated via cal¬ 
ibration measurements. The fraction F = 0.52 ± 0.03 
is constant over the considered energy range. The cti 
for each ROI is calculated from the resolution function 
shown in [Figure 1| The values range from 7.2 keV up to 
12.8 keV, and the corresponding uncertainties range from 
0.09 keV to 0.15 keV. 

All efficiencies are estimated as the ratio of events 
which fulfill the corresponding selection criteria and the 
total number of events. The depth requirement efficiency 
is evaluated with a pulse shape simulation. It is based 
on a simulation of the electric fields for the coplanar elec¬ 
trode layout which simulates charge propagation includ¬ 
ing diffusion and electromagnetic repulsion of the charge 
cloud. To estimate possible systematic deviations from 
the data, the results of the simulation are compared to 
an efficiency calculated using the observed ^^^Cd spec¬ 
trum with an endpoint of 320 keV. For further informa¬ 
tion about the spectrum see Ref. [16]. ^^^Cd is homoge¬ 
nously distributed inside the detector and gives a good 
measure of the efficiency losses caused by the depth re¬ 
quirement in an energy range from 180 keV to 300 keV. 
The absolute difference between both efficiency measure¬ 
ments is about 0.1. A possible source of these systematic 
deviations are distortions of the energy reconstruction 
near the anode grids. The efficiency is calculated for 
each ROI separately and range from 0.74 (^^^Cd) to 0.78 
(^^®Cd). The associated absolute systematic uncertainty 
is ±0.1. 

The impact of the LSE requirement on the efficiency 
is estimated with the double escape peak of the 
7 line of the calibration measurements. This is advan¬ 
tageous because the early rise time also depends on the 
spatial distribution of the energy deposits. This distri¬ 
bution is smaller for the double escape peak and QvjSfi 
decay events and it is more confined than events caused 
by multiple-scattered photons. The efficiency is system¬ 
atically underestimated even if calculated at the double 
escape peak because of the underlying background caused 
by The emitted gamma line overlaps strongly with 

the double escape peak and hence the fraction of the 
^^^Bi line is minimized by using the low energy side of 
the peak. Based on measurements with calibration data 
the efficiency is estimated to be 0.76 ± 0.02. 

The intrinsic efficiency is the probability to deposit the 
whole energy of a Qvj3j3 event inside the active volume of 
one detector and is estimated with Monte-Carlo simula¬ 
tions. The values range from 0.62 to 0.96. The abso¬ 
lute uncertainty of the calculated efficiencies is not ex¬ 
pected to be larger than 0.002. The simulation is based 
on GEANT4 9.6, which models the particle propagation 
and interaction with matter |23| . The energy and an¬ 
gular distributions of the physics processes are predicted 


by an event generator based on DECAYO |2Tj. To es¬ 
timate detector specific efficiency losses, causes, e.g., by 
dead layers, each detector is scanned with a ^^^Cs point 
source with a known activity in three different distances. 
The peak count rates of these three measurements are 
fitted with f{r) = a/r^. The fit parameter a is pro¬ 
portional to the activity of the source, the solid angle 
fraction covered by the detector, the emission probabil¬ 
ity of the 7 -line and the efficiency of the detector. For 
the estimation of the efficiency of the full detector array, 
the detector with the highest efficiency is assumed to be 
fully efficient. All other detector efficiencies are calcu¬ 
lated in relation to this detector. The average detection 
efficiency of the detector array is 0.89 ± 0.01. 

TABLE I: The first part lists all means and widths 
which are used as prior informations. The second part 
lists all values which are used to calculate the total 
uncertainties. 





i 28 Te 

"°Zn 

isoTe 


Q-value 

[keV] 

543 

867 

997 2528 

2814 

O' Q-tot 

[keV] 

2.1 

2.1 

2.8 

1.9 

2 

F 


0.52 

0.52 

0.52 

0.52 

0.52 

ap 


0.03 

0.03 

0.03 

0.03 

0.03 

Ol 

[keV] 

7.23 

8.27 

8.65 

12.29 

12.86 

0 1 ,unc 

[keV] 

0.1 

0.09 

0.09 

0.15 

0.15 

£tot 


0.48 

0.46 

0.45 

0.34 

0.33 

^etot 


0.07 

0.07 

0.07 

0.05 

0.05 

N 

[1q 23 atoms j 

6.59 

8.08 

0.015 

8.62 

1.73 

(JN 

^1q 23 atoms j 

0.132 

0 

0.003 

0 

0.035 

O Q-value 

[keV] 

0.9 

0.9 

2.1 

0.01 

0.13 

O cal 

[keV] 

1.9 

1.9 

1.9 

1.9 

1.9 

^depth 


0.74 

0.74 

0.74 

0.76 

0.78 

'^Edepth 


0.1 

0.1 

0.1 

0.1 

0.1 

^Ise 


0.76 

0.76 

0.76 

0.76 

0.76 



0.02 

0.02 

0.02 

0.02 

0.02 

^int 


0.96 

0.92 

0.90 

0.66 

0.62 



0.002 

0.002 

0.002 

0.002 

0.002 

^det 


0.89 

0.89 

0.89 

0.89 

0.89 



0.01 

0.01 

0.01 

0.01 

0.01 


The total efficiency uncertainty is obtained by error 
propagation and is shown in [Table lllj Runs with the 
highest trigger threshold of 550 keV are excluded for the 
^^"‘Cd analysis and are well below all other region of in¬ 
terests. Therefore the trigger efficiency is assumed to be 
100 %. 

One virtual CdZnTe molecule consists of 
Cdo.gZng.iTei. This composition is consistent with 
inductively coupled plasma mass spectroscopy mea¬ 
surements, which were done with a CdZnTe sample. 
However, the fraction of Zn can vary throughout the 
crystals and therefore a systematic uncertainty of 
±2% points is introduced for the Cd fraction, which 
corresponds to a variation of ±20 % in the Zn fraction. 

The detector lifetime is measured by the internal 
FADC clock. The uncertainty of each clock is of the 
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order of 10^. The mass of each detector is measured 
before installation with an accuracy of milligrams which 
is negligible compared to a single detector mass of approx 
6g. Therefore, the exposure uncertainty is neglected in 
this analysis. 

IX. RESULTS 


TABLE III: Results of the Bayesian signal estimation 
after folding all systematic uncertainties. The second 
column shows the natural abundance of atoms per kg 
which is used for the signal estimation, for its 
uncertainty see |TableT| In the third column the 
background index b for the different ROIs is presented. 
The fourth column reports the limit at 90 % credibility 
and the last column reports the calculated Bayes factor. 


The energy spectrum is fitted in several ROIs for the 
isotopes listed in [Table Ilj The table lists the exposure 
after the run selection, the region of interest, the Q-values 
|22| and the interval of the uniform prior used for each 
isotope. The energy spectra including the best fit model 
and the corresponding uncertainty bands are shown in 
Figs. |4a|to|4^ 

TABLE II: Regions of interests and the corresponding 
Q-values and prior probabilities for each isotope. The 
priors vary because of different isotopic abundancies 
and efficiencies 


Isotope 

A 

[kgd] 

ROI 

[keV] 

Q 

[keV] 

^1/2 prior 
[yr-'l 


212.8 

483 to 584 

543 

0 to 2 x lO-rs 

i28Te 

216.1 

809 to 926 

867 

0 to 2 x 10 “^® 

"°Zn 

216.1 

940 to 1062 

997 

0 to 1 x lO-re 

i30Te 

216.1 

2440 to 2671 

2528 

0 to 4 x 10“®° 


216.1 

2723 to 2904 

2814 

0 to 2 x 10 “r° 


For each ROI, a signal model. Hi, and a background 
only model, Hq, are defined. The Bayes factor is defined 
as 


P{D\Hoy 


(7) 


Isotope 

A/10®° 

b 

Ti/2 90% C.L. 

K 


[ atoms 1 

1 kK J 

\ counts I 

1 keV kg yr J 

[10®Vr] 


ri^Cd 

6.59 

213.9t; 7 

1.6 

0.07 

i28Te 

8.08 

65.5t°:® 

1.9 

0.17 

roZn 

0.015 

45.1+?° 

6.8 X 10“° 

0.06 

laoTe 

8.62 

0 ft+O-l 
'^•0-0.3 

6.1 

0.14 

rr°Cd 

1.73 

9 y+0.1 
' -0.2 

1.1 

0.27 


The fit of ^°Zn prefers zero events which results in a limit 
of Ti /2 > 6.8 X 10^® yr. The included in the 

^®*^Te fit yields zero events, so there is no evidence for any 
background contribution from the Thorium decay chain. 
The limit for ^®°Te is T 1/2 > 6.1 x 10^^ yr. In the ^^®Cd 
ROI the fit yields zero events for the 0^/3/3 decay mode. 
The resulting limit is T 1/2 > IT x 10^^ yr. Switching to 
a prior uniform in the effective Majorana neutrino mass 
increases the limits on the half-life by 30-40%, depending 
on the isotope. 

[Table III] reports the abundance per kg detector mass 
of the analyzed isotope, the total efficiency, the back¬ 
ground level in each ROI and the 90% credibility lower 
limit on the half-life for each isotope. 


X. CONCLUSION 


where P{D\Hi) is the probability of the data given 
the signal model and P{D\Hq) the probability of the 
data given the Hq model. Bayes factors below 1 indi¬ 
cates that the signal hypothesis is disfavored against the 
background-only hypothesis. 

All Bayes factors are well below 1, and thus the 
signal hypothesis is disfavored for each isotope. The 
absolute values of K are listed in ITable lITl The 
observed events in the ROIs range from 288 (^^®Cd) 
to 18 670 (^^^Cd). The corresponding background in¬ 
dices range from 2.7counts/(keVkgyr) (^^®Cd) up to 
213.9counts/(keVkgyr) (^^"‘Cd). 

For ^®°Te the 2°®T1 line at 2 614keV is included in the 
analysis to scan for possible contributions of the thorium 
decay chain. If the signal hypothesis is rejected, a lower 
limit on the half life at 90 % credibility is calculated. For 
the analysis results in a lower limit of T 1/2 > 
1.6 X 10^^ yr. The included annihilation line has a signal 
strength of 1770 counts. For ^^®Te the analysis yields 
zero events resulting in a limit of T 1/2 > 1-9 x 10^^ yr. 


A search for Oiz/3/3 decay with the COBRA demonstra¬ 
tor setup has been performed using data corresponding to 
an exposure of 234.7 kg d. No indication for 0vPl3 decay 
have been found and 90 % credibility limits have been set 
for different isotopes. The COBRA demonstrator contin¬ 
ues to collect data. A layer of large segmented CdZnTe 
detectors will be installed in the near future, which is 
expected to increase the sensitivity considerably. 
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